Globular Clusters in Fornax: Does M° Depend on Environment? 

John P. Blakeslee 1 and John L. Tonry 

Department of Physics, 6-216, Massachusetts Institute of Technology, 77 Massachusetts Ave., 

Cambridge, MA 02139 



ABSTRACT 

We present the y-band globular cluster luminosity functions (GCLFs) of the Fornax 
Cluster galaxies NGC 1344, NGC 1380, NGC 1399, and NGC 1404. Our observations 
reach to V = 24.8, roughly one magnitude beyond the GCLF turnover m v , with ~ 90% 
completeness. From the amplitude of the galaxy surface brightness fluctuations, we also 
estimate the number of globular clusters fainter than this cutoff magnitude. The GCLFs 
of these galaxies are well fitted by Gaussians; the weighted means of their turnover 
magnitudes and dispersions are (m v ) = 23.88 ± 0.10 mag and (a) = 1.35 ± 0.07 mag. 
The assumption of a universal value for the absolute magnitude of the turnover M v 
places the Fornax cluster 0.13 ± 0.11 mag more distant than Virgo. However, in light 
of recent Cepheid and other high-precision distance measurements, as well as ongoing 
HST observations of GCLFs for the purpose of determining the extra-galactic distance 
scale, we choose to re-examine the universal GCLF hypothesis. Based on data from 
groups and clusters of galaxies, we find evidence that M v becomes fainter as the local 
density of galaxies increases. We speculate on the possible cause of this trend; if it is 
confirmed, GCLF observations will be less useful for determining distances, but may 
provide important information for constraining theories of star formation in primordial 
galaxy halos. 

Subject headings: galaxies: clusters: individual (Fornax) — galaxies: distances and 
redshifts — galaxies: star clusters — globular clusters: general 
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1. Introduction 

The globular cluster luminosity function 
(GCLF) is often employed standard-candle 
distance indicator based on the assumption a uni- 
versal value for its mean, or turnover, magni- 
tude M° (see |Jacoby et al. 1992| for a review 



of the method). Until recently, it was impossi- 
ble to apply the GCLF method to determine the 
distances of galaxies further away than Virgo. 
Now, with HST and improvements in ground- 
based seeing and instrumentation, it becomes po- 
tentially much more powerful for determining the 
extra-galactic distance scale. Furthermore, new 
Cepheid and other high-precision local distance 
measures would allow for a firm calibration of the 
method. 

Baum et al. (1995a,b) have used HST to ob- 
serve the the globular clusters (GCs) of the 
Coma galaxies NGC 4881 and IC 4051 down to 
V = 27.6 and V = 28.4, respectively. They de- 
rive values of the Hubble constant Hq near 60 
km/s/Mpc. Also with HST, Whitmore et al. 
(1995) studied the GCs of the extremely rich M87 
system to two magnitudes beyond My (the V- 
band GCLF turnover) and derived Hq = 78 ± 11 
km/s/Mpc. Only a small part of the discrep- 
ancy in derived Hq values can be accounted for 
by the different calibrations used by the two 
groups. This situation leads one to suspect that 
the GCLFs themselves may be intrinsically dif- 
ferent, especially as there remains no firmly es- 
tablished physical basis for assuming a univer- 
sal M . Previously, there have been suspicions 
that M° is different for spirals and ellipticals 
((Seeker fc Harris 1993 ), with the root cause of 
this difference being metallicity variations (Ash- 
man, Conti, & Zepf 1995), but M v was assumed 
not to vary among large ellipticals. As Whitmore 
et al. candidly remark, this "crucial assumption" 
of a universal GCLF is "a hypothesis that needs 
further verification." 

In this Letter, we examine the current state 
of the universal GCLF hypothesis. First, we 



present new observations of GCs around four 
Fornax galaxies. Fornax is an important clus- 
ter for testing distance determination methods, 
as it is spatially much more concentrated than 
Virgo while being at nearly the same distance 
(e.g. [Tonry 1991 ; Ciardullo, Jacoby, & Tonry 
1993). We find that the GCLF exhibits remark- 
ably little variation for these Fornax galaxies. 
Next, we use independent distance measurements 
to galaxies and galaxy groups to compare derived 
My values in different environments. We find 
somewhat startling evidence that My becomes 
fainter as the local density of galaxies increases. 
Further verification is once again needed, but if 
the observed trend proves real, it would have ma- 
jor implications for the GCLF method of distance 
measurement as well as for theories of GC for- 
mation. We conclude with a discussion of these 
implications, in particular how the local galaxy 
environment may govern the properties of GC 
populations. 

2. Observations and Reductions 

We observed the Fornax Cluster galaxies NGC 
1316, NGC 1344, NGC 1380, NGC 1399, and 
NGC 1404 in 1995 August with the Tek 2048 2 #4 
CCD detector at the Cassegrain focus of the 
4 m telescope at Cerro Tololo. Four 600 s V- 
band exposures were taken of each galaxy, ex- 
cept NGC 1316, for which five 600 s exposures 
were taken. We also obtained 2400 s of integra- 
tion on a background field 1?5 west of the cD 
NGC 1399. The image scale was Of! 158 pix" 1 ; 
however, the chip was slightly vignetted around 
the edges, and we shifted the telescope ~ 6" 
between individual exposures, so the final field 
size which received the full integration time was 
about 5'.1 x 5fl. We processed the images as de- 
scribed by Ajhar, Blakeslee, & Tonry (1994) and 
Blakeslee & Tonry (1995; hereafter BT95). The 
seeing in the final NGC 1399 image was / .'94, 
while the seeing in the other images ranged from 
l'/03 to 1'.'05. The photometry was calibrated us- 
ing Landolt (1992) standard stars; there is no 
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detectable Galactic extinction in the directions 



of these galaxies flBurstein k, Heiles 1984) . 

After subtracting smooth models of the galaxy 
surface brightness profiles, we used a version 
of the program DoPHOT (Schechter, Mateo, 
& Saha 1993) for the point source photome- 
try. Completeness corrections were determined 
by scaling and adding grids (so as to avoid ar- 
tificial crowding) of 32 x 32 "cloned" PSF stars 
and then finding them again with DoPHOT. The 
scaling was done at 0.4 mag intervals; we inter- 
polated to find the completeness corrections at 
intermediate magnitudes. The images were then 
divided up into three radial regions: 20-40", 40- 
80", and 80-160". We settled on a cutoff mag- 
nitude m c of V = 24.4 for the innermost region 
in each galaxy; at this magnitude, the complete- 
ness levels for this region ranged from 77% to to 
88%. For the intermediate region, the complete- 
ness levels were in the same range at V = 24.8, 
so we used this for the cutoff magnitude. For 
simplicity, we also used V = 24.8 for m c in the 
outermost region, although the completeness lev- 
els for this region ranged from 90% to 96% at 
this magnitude. The photometric error is typ- 
ically ~ 0.12 mag at V = 24.8. Any objects 
classified by DoPHOT as extended were excluded 
from further analysis. Objects brighter than m c 
in each region were then binned in magnitude 
and our completeness corrections were applied. 
We subtracted the completeness-corrected lumi- 
nosity function of the unresolved objects in the 
background field from the luminosity functions 
of the objects in the program fields to produce 
the final GCLFs presented below. 

After removing all objects brighter than m c 
in each region, we measured the PSF-convolved 
variance, or fluctuations, remaining in the image. 
Our variance analysis method is described in de- 
tail by Tonry et al. (1990). The variance mea- 
surement effectively acts as an "extra bin" in fit- 
ting the GCLF. BT95 demonstrate how this mea- 
surement can be used for deriving GCLF param- 
eters. The conversion from measured variances 



to GC densities is done in the same way here, 
but the GC counts are treated differently in that 
they are binned for more GCLF shape informa- 
tion. In addition, we leave free param- 
eter, instead of varying it only within some pre- 
supposed acceptable range, as in BT95. A more 
detailed account of our point source photome- 
try and positions, completeness experiments, and 
fluctuation measurements will be provided else- 
where ( [Blakeslee et al. 1996j ). 

3. Results 

We defer the analysis of the spatial structure 
and total sizes of the GC populations and concen- 
trate on the luminosity functions. We also defer 
any further discussion of our observations of GCs 
in the giant disturbed galaxy NGC 1316 (Fornax 
A). The GCLF of this galaxy was not well fitted 
by a Gaussian model (x 2 ^> 3 instead of ~ 1) and 
would require a more thorough analysis than we 
can provide here. 

3.1. The GCLF in Fornax 

Figure 1 presents the l/-band GCLFs of the 
four Fornax galaxies. The plotted curves are 
Gaussians having a and my values which are the 
weighted means of the values found in our three 
analysis regions of each galaxy. The values were 
derived in the separate regions by x 2 minimiza- 
tions using the counts brighter than the individ- 
ual cutoff magnitudes and our variance measure- 
ments. For purposes of displaying the GCLFs, 
however, we used m c = 24.8 everywhere, applied 
our completeness corrections, then binned the re- 
gions all together. 

Table 1 lists our final values for the GCLF 
parameters; they are nearly identical within the 
errors. The cD NGC 1399 is the only one with a 
well-studied GCLF. Geisler & Forte (1990) found 



V 



23.45 for this galaxy, but assumed a 



1.20 mag (these parameters are correlated when 
the limiting magnitude is near the turnover). 
Bridges, Hanes, <fe Harris (1991) used a ~ 1.40, 
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Table 1 

Final Fornax Gaussian GCLF parameters. 



Galaxy 


Type 


V T 


a 


± 




± 


N 


± 


N1344 


E5 


10.35 


1.35 


0.18 


23.80 


0.25 


140 


20 


N1380 


SO 


9.91 


1.30 


0.17 


24.05 


0.25 


375 


35 


N1399 


El/cD 


9.57 


1.38 


0.09 


23.83 


0.15 


1120 


45 


N1404 


El 


9.98 


1.32 


0.14 


23.92 


0.20 


300 


30 



Note. — Columns list: galaxy name, Hubble type, total apparent V-band magnitude ( de Vaucouleurs et al 
19E[l]), Gaussian dispersion of the GCLF, ]/-band GCLF turnover magnitude, and total number of GCs which went 
into the Gaussian fits (following corrections for incompleteness and background). 
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The glubulai clusLei liiniinusiiv imiciiuiis 

of the four Fornax galaxies. Filled symbols represent 
the final counts following incompleteness and back- 
ground corrections. Open symbols show what the 
counts would be with background subtraction, but 
no incompleteness corrections. Errorbars represent 
the uncertainties in our corrections and the Poisson 
errors in the counts. 



and found m v = 23.85 ±0.30, in close agreement 
with our value. 

As weighted means of the Fornax GCLF pa- 
rameters, we take {a) For = 1-35 ± 0.07 mag; 

23.88 ±0.10. Whitmoreet al. (1995) 
= 23.72 ± 0.06 for M87. We aver- 
age this with the values for NGC 4472 and NGC 
4649 QSecker fc Harris 1993|) using (B — V) = 
0.75 ± 0.05 QCouture et al. 199q , 1991) to find 
(my)yi r = 23.75 ± 0.05. Thus, assuming a uni- 
versal M v yields A(m - M) = 0.13 ± 0.11 for 
the relative Fornax- Virgo distance modulus. If 



("iy) For = 

found my 



we were to include N GC 4365 (pecker fc Harris 
1993|) and NGC 4636 flKissler et al. 1994| ) in the 



Virgo average, the relative modulus would drop 
by ~ 0.1, but both these galaxies are question- 



able Virgo members (Tully 1987b; Tonry et al 



1990) 



3.2. Does M° Depend on Environment? 

We now forsake the assumption of universal- 
ity and compare the Fornax GCLF with those 
observed elsewhere. To do this, we need a self- 
consistent set of independent distance determina- 
tions to galaxies or groups in which m v has been 
measured. We start by fixing (m — M) = 31.0 
as the Virgo distance modulus; this is both the 
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Jacoby et al. (1992) value and the latest HST 



Cepheid Key Project result flFreedman et al 
1991]). For the relative Fornax- Virgo distance 
modulus, we take A(m — M) = 0.25 ± 0.08 from 
an average of the PNLF, SBF, D n - a, SN la, 
and Tully-Fisher methods (Ciardullo et al. 1993; 
Faber et al. 19891 ; |Riess 1996fc |Willick 19961 )" 
We use measurements of m v for M31 ( pecker 



19$|, M81 ( Perehnuter fc Racine 1995) ), the Leo 
group ellipticals ( Harris 1990 ), and the HST limit 



on my for NGC 4881 in Coma ( Baum et al 



19£5a). The Cepheid distance moduli to M31 
and M81 are 24.43 ± 0.10 and 27.80 ± 0.20, re- 



spectively ( Freedman &: Madore 199C ; Freedman 
et 3,1. 1994 ). This M31 distance modulus is the 
proper one to use, as it is the one assumed by 
Jacoby et al. (1992) and is consistent with the 
more recent HST Cepheid distances. For the 
Leo group, we average the HST Cepheid distance 
to the spiral M96 ( Tanvir et al. 1995| ) with the 
PNLF and SBF distances to the ellipticals ( |Ciar- 
du jlo et al. 1993j ) and get (m-M) =30. 2 ±0.13. 



The relative Virgo-Coma distance modulus is 
well determined at A(m — M) = 3.71 ±0.10 (e.g. 
van den Bergh 1992^ Whitmore et al. 1995). Fi- 
nally, in an effort to preserve neutrality in the 
controversy over the RR Lyrae calibration (see 
van den Bergh 1995| ), we omit the Milky Way 
from our discussion, remarking only that recent 
MW M v values have ranged from -7.29 ± 0.13 
( pecker 1992Q to -7.60 ±0.11 ( |Sandage fc TkT^- 
ma |nn 19*95 ; see also the discussion by Baum et 
al " 
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Figure 2 shows the resulting My values plot- 
ted against the the velocity dispersions of the 
groups and clusters, from Tully (1987a) and 
Zabludoff, Huchra, & Geller (1990). We use ve- 
locity dispersion as the most convenient indica- 
tor of the depth of the local potential well; it 
closely correlates with Tully 's estimated group 
densities and with cluster richness. There is a 
trend of decreasing turnover luminosity with in- 
creasing local density. The offset in My between 
the small groups and Fornax/Virgo is 0.4 mag. 



The use of the straight Cepheid distance to Leo 
would move its My brighter by 0.1 mag, further 
away from the other ellipticals; the inclusion of 
NGC 4365 and NGC 4636 would move the Virgo 
My fainter. In addition, the preliminary result 
m v w 28.0 {My rj -6.7) for IC 4051 (Baum et 
al. 1995b) indicates that there may be another 
~ 0.5 mag offset in My for the very rich Coma 
cluster. Thus, we believe we are seeing real evi- 
dence for an environmental dependence of My. 

4. Discussion 

We have found that the GCLF is remarkably 
constant within the Fornax cluster, but, as Fig- 
ure 2 shows, varies with environment. Ashman 
et al. (1995) suggested that metallicity differ- 
ences result in My values which are systemat- 
ically brighter by ~ 0.15 mag for spirals. Large 
ellipticals usually do have higher metallicity GC 
populations; however, NGC 4881 in Coma has a 
GC color/metallicity distribution similar to that 
of the MW (|Baum et al. 1995a|) , yet its My is 
very faint. In addition, the Leo elliptical NGC 
3379, with its relatively high metallicity GC pop- 



ulation (Ajhar et al. 1994), has an exceedingly 
bright My , though with a large uncertainty ( Har- 



ris l99C| ). Finally, we note that the magnitude of 



the environmental effect we propose is a factor of 
3-6 larger than the Ashman et al. My metallicity 
shift. 

We suggest that the most straightforward way 
to produce the present-day near-Gaussian GCLF 
is to assume that two simultaneous and compet- 
ing effects were operating when GCs formed: a 
"creation" process which preferentially created 
low-mass GCs, cutting the mass function off at 
the high end, and a "destruction" process which 
inhibited the formation of, or quickly destroyed, 
low-mass GCs. If each process operated in a 
manner which was independent of the details of 
the environment, then the final mass (luminos- 
ity) function would be universal, but if one de- 
pended more sensitively on environment than did 
the other, the final mass function would vary. 
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Fig. 2.— The GCLF turnover magnitude My plot- 
ted against the velocity dispersion of the host galaxy's 
environment, used as a measure of the local density. 
Sec text for details. 
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Fig. 3. — The effects of variable GC creation/de- 
struction mechanisms. In part (a) the GC creation 
process, shown as a power-law growing to smaller 
mass, is universal (dark solid line), and the destruc- 
tion process, a power law which wipes out low-mass 
objects, varies with environment (dashed lines). Since 
the total number of GCs is the integral under the 
intersecting creation/destruction lines, this situation 
results in an anti-correlation between GC specific 
frequency Sn and GC mean logarithmic luminosity 
— M°. In part (b) the creation process (dashed lines) 
varies with environment, while the destruction pro- 
cess (dark solid line) is universal. The result here is 
a positive correlation between Sn and — M°. 



This situation is schematically illustrated in Fig- 
ure 3, where we use Sn for the GC "specific fre- 
quency" (number of GCs per unit luminosity of 
the host galaxy). In Figure 3a, we assume that 
the creation process is relatively universal, but 
that the destruction process varies. This leads 
to M° being variable, and predicts an inverse 
correlation between the number of GCs formed 
and their mean brightness. On the other hand, if 
the destruction process is constant and the cre- 
ation process is more variable we will again get 
a variable M°, but with a direct correlation be- 
tween the number of GCs and their luminosity, 
as shown in Figure 3b. 

Empirically, we think we see evidence for the 
latter sort of behavior among "coeval" galaxies, 
i.e., those located within the same physical as- 
sociation. In Virgo, for instance, M87 has a 
very large Sn and a slightly brighter M° than 
its close neighbors, and similarly in Fornax for 
NGC 1399. In this context, the GCLF would 
depend on the extent to which the host galaxy 
dominated its local environment. On the other 
hand, the main point of this paper is that we see 
the former behavior among very "heterogeneous" 
systems of galaxies. Young groups dominated by 
spirals have fewer GCs than galaxy clusters such 
as Virgo, which in turn may have fewer than rich 
clusters such as Coma, and we find that the cen- 
tral luminosity of the GCLF is declining along 
this sequence. 

As an example of how such an interplay of op- 
posing processes might work in practice, we con- 
sider the common picture of structure formation 
through gravitational instability. Here, the "cre- 
ation" process is the primordial spectrum of den- 
sity fluctuations which favors low-mass clusters, 
and the "destruction" process is the inhibition of 
the collapse of low-mass objects resulting from 
the Jeans mass. In this picture, the Jeans mass 
can be a very rapidly growing function of time 
(Tegmark et al. 1996), and the densest systems 
of galaxies, forming first, would have experienced 
a less restrictive low-mass cutoff and hence have 
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more, and fainter, GCs. This is precisely the case 
depicted in Figure 3a. 

Harris k Pudritz (1994) have proposed a de- 
tailed astrophysical theory of GC formation which 
is perhaps more illustrative of the "coeval" case 
of Figure 3b. They suggest that the "creation" 
process is made more efficient by the larger ex- 
ternal pressures of dense environments. Their 
"destruction" process is the tidal disruption and 
evaporation of low-mass GCs, and this might be 
less sensitive to environment (although see Mu- 
rali k Weinberg 1996). Of course, if the cutoff 
is very abrupt (a steep "destruction" line), then 
M° will not correlate very strongly with Sjy. 

We do not have better "creation" and "de- 
struction" processes to offer than have been ad- 
vanced elsewhere, but we believe that this de- 
scription is a profitable way to frame the dis- 
cussion. Until now, the assumed constancy of 
M° has been a serious obstacle to reasonable 
models for GC formation, so we conclude by re- 
emphasizing our primary point. The GCLF ap- 
parently does depend on environment, with My 
being fainter in denser regions, although it may 
be remarkably constant within a single group of 
galaxies. This dependence will present challenges 
for the use of the GCLF as a distance indicator. 
On the other hand, it opens the door for correla- 
tions between M° and Sn, and M° and environ- 
ment, which may yield valuable insights into the 
conditions and processes which prevailed at the 
time of GC/galaxy formation. 

J.B. thanks Ed Ajhar for invaluable help while 
observing, as well as profitable and sometimes 
overly animated discussions, and the CTIO 
support staff for their tireless troubleshooting 
efforts. We also gratefully acknowledge helpful 
discussions with Lam Hui and Paul Schechter. 
This research was supported by NSF grant AST94- 
01519. 
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